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Abstract — We examined dipolar field effects on the microwave 
generation in the domain wall spin-valve by solving simultane- 
ously the Landau-Lifshitz-Gilbert and Zhang-Levy-Fert diffusion 
equations. By numerically analyzing dipolar field dependence, 
we showed that the microwave generation needs the dipole- 
dipole interaction for a 180° domain wall and the amplitude of 
the microwave voltage signal depends strongly on the exchange 
length. In order to design a microwave generator using the 
domain wall spin-valve with high efficiency, we propose that the 
materials with short exchange length are preferred. 

Index Terms — CPP-GMR, current-confined-path, micromag- 
netic simulation, domain wall 

IN recent years, much effort has been devoted for devel- 
opment of a microwave generator using a nano-scale spin- 
transfer torque oscillator [Q, ED, El, El, flU, (6), Q, (8). As a 
powerful candidate of the nano-scale spin-ttansfer oscillator, 
a domain wall spin- valve is proposed 0. The domain wall 
spin-valve consists of a magnetic nano-multilayer, the main 
part of which is a nano-oxide-layer (NOL) sandwiched by 
ferromagnetic layers iflOl . The NOL contains a lot of fine 
contacts filled up with ferromagnetic metal. The contacts 
enable dc-current to be conducted by bias voltage. Under the 
dc-current, a microwave voltage signal is switched on when 
domain walls in the contacts are constructed. This implies that 
the microwave voltage signal is considered to originate from 
dynamics of the domain walls driven by the dc-current in the 
contact. 

The model of the domain wall spin valve is shown in 
Fig. 0I a )- Based on the model we previously reproduces the 
microwave generation in the domain wall spin-valve under dc- 
current by a micromagnetic simulation technique ifTTl . In the 
simulation the spin-transfer torque drives a uniform rotation 
of local magnetizations as shown in Fig. Hlb). In the rotation 
the magnetic structure oscillates between the Neel and Bloch 
walls. The oscillation generates the microwave voltage signal 
by inducing an oscillation of spin accumulation which con- 
tributes voltage drop. However the reason why the oscillation 
of the magnetic structure induces the microwave voltage signal 
is still unclear because the voltage drop induced by the spin 
accumulation mainly depends not on the magnetic structures 
but on the thickness of the domain wall. The mechanism 
should be clarified in order to design a microwave generator 
using the domain wall spin valve with high efficiency. 

In the present work, we propose a scenario of a breathing 
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Fig. 1. The schematic picture of the model of a geometrically confined 
domain wall (a) and the rotation of the magnetic structure (b). White and 
black arrows denote the magnetic structures for the Neel and Bloch walls, 
respectively. The dotted circle denote the rotation of the local magnetization 
configuration. 



mode excitation in the domain wall, which relates the rotation 
of the magnetic structure with the oscillation of the spin accu- 
mulation and is induced by both the dipole-dipole interaction 
and spin-transfer- torque. The scenario is as follows: The spin- 
transfer torque induces the rotation of the magnetic structure 
IfTTl . ifPHl . Ifl3ll . Then a breathing mode, which is an oscillation 
of the thickness of the domain wall, appears because the 
dipole-dipole interaction effectively introduces a difference 
between the thickness of the Neel and Bloch walls. As a result, 
the dc-current is converted to an ac voltage signal because the 
resistance of the domain wall due to the spin accumulation is 
proportional to squared inverse of its thickness |[l4l . 

The scenario predicts strong dependence of magnitude of 
the microwave voltage signal on the exchange length, l cx , 
defined by yj Jdd/2-fr'd, where Jdd and are coupling con- 
stants of the exchange and dipole-dipole interaction between 
local magnetizations, respectively, because the dipolar field 
yields the breathing mode. In the present paper, to confirm 
the prediction of our scenario, we simulated dynamics of 
the local magnetization configuration and spin accumulation 
in the domain wall spin valve by solving simultaneously 
the Landau-Lifshitz-Gilbert lfT31l and Zhang-Levy-Fert diffu- 
sion equations! 16i| an d investigated the dipolar field effects by 
artificially controlling the exchange length. We considered the 
values of the exchange length of CoFe, Permalloy and infinity. 
Based on numerically analyzing exchange length dependence, 
we showed that the efficient microwave generation needs the 
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strong dipole-dipole interaction for the 180° domain wall 
geometrically confined in the contact fTTI and the magnitude 
of the microwave voltage signal depends strongly on the 
exchange length. Our results suggests that the materials with 
short exchange length is preferred in order to develop the 
microwave generator using the domain wall spin-valve with 
high efficiency. 

The system consists of two magnetic electrode layers and a 
non-magnetic insulator layer sandwiched by the electrodes as 
shown in Fig.[TJa). The insulator layer contains magnetic con- 
tact which geometrically confines a domain wall. The system is 
divided into about 1500 hexahedral finite elements, where spin 
accumulation and demagnetization field are evaluated. The size 
of magnetic metal electrode at the bottom and top layers is 8 
x 8 x 4 nm. The shape of the contact between two electrodes 
is a rotated elliptic arch around the center axis perpendicular 
to the layers, the diameters of the bottom and center layers 
of which are set at 6 and 4 nm, respectively. The size of the 
system is larger than the diffusion length of 2nm of CoFe lfT8l 
at room temperatures. Thus the electronic system is diffusive 
at least. We deal with the electron system in diffusive limit, 
unlike a previous ballistic treatment! 19). 

The local magnetizations are expressed by the classical 
spins on the simple cubic lattice with the lattice constant of 
a = 0.4nm. The Hamiltonian TL is given by 

H = - Jdd 22 Si ■ Sj + J s d ^ Si ■ brn l 



(id) 



-f*I>-/ d H |^|3 3^ 



S(r), (1) 



where fi represents the relative coordinate of the i-th site 
from the position r; Si, the classical Heisenberg spin with 
absolute value of unity and 5rhi, local spin accumulation 
density at the z-th site. The first term in the right hand 
side of Eq. (Q]) expresses the exchange interactions between 
local magnetizations at nearest neighbor sites. The exchange 
coupling constant denoted by Jdd is related to the exchange 
stiffness constant in continuous limit, A, and a lattice constant, 
a, with Jdd = 2aA. We fix Jdd at 0.04 eV which is of the 
order of the transition temperature of the typical material, 
T c ~ 0.8(Permalloy)-1.2(CoFe) xl0 3 K, and the material 
dependence of Jdd is neglected because the dependence is 
much weaker than other material parameters and because the 
spin-transfer torque depends only on Jdd through the small 
correction of the order of the small parameter defined by 
11 JddT«f ll20l . where r s f is a spin relaxation time by a spin- 
orbit coupling. 

The second term expresses the s-d exchange interactions 
between local magnetizations and the spin accumulation at 
each site ll20l . The s-d exchange coupling constant J s d is set 
at 0.1 eV in accordance with Ref. |[T6l . The term reproduces 
the spin-transfer torque induced by the spin accumulation. In 
the present paper, we adopted the Zhang-Levy-Fert diffusion 
equations |[T6l for the evaluation of the spin accumulation. 

The third term denotes the dipole-dipole interaction energy. 
In order to deal with the dipole-dipole interaction in a complex 
system shape, we adopt a finite element - boundary element 



(FEM-BEM) hybrid method ED on the spin field on the 
continuum space, S(r), which is defined for each position, 
r, by the spin field extrapolated from the lattice sites to the 
position. The exchange length for the dipole-dipole interaction 
is l ex = \J Jdd/2Kd nm. We examined the cases of the 
exchange length of 3nm(CoFe), 5nm(Permalloy) and 00 by 
controlling artificially Kd in order to investigate dipolar field 
effects. 

In this paper we assume the system size of 4nm is larger 
than the mean free path. In this case, the local spin accumu- 
lation density. Sfhi, is determined by solving the following 
Zhang-Levy-Fert diffusion equation] 16|, 



— 5rh(f 



V {,8S(r)l(r) + A(S(r))Sm(f)} 
+ J f ±Sm(r)xS(r)+ A ™^, (2) 
A(S(f)) = 2D \l - p 2 S(r) ® S(r)v] , (3) 



where j e , Co, Dq, (3 and r denote electronic current, con- 
ductivity, diffusion constant, polarization of resistivity and 
relaxation time due to a spin-orbit interaction respectively. 
Equations (fSJl-Ol) are solved numerically with combining con- 
tinuous equation for electronic current, Vj e = 0. Co and 
(3 are taken to be those for conventional ferromagnets as 
Co =70£!nm and j3 = 0.65, respectively. Do and r are ob- 
tained, respectively, by the Einstein relation, Co = 2e 2 NpD 
and A = ^/2r£>o(l — P 2 ) for given spin diffusion length A, 
density of states at the Fermi level, N-p and electron charge e. 
Here we employ A=12nm and N-p =7.5nm~ 3 eV~ 1 . 

We adopted an adiabatic approximation where the relaxation 
time of the spin accumulation is much shorter than that 
of the local magnetization. In that case, we can neglect 
dynamics of the local magnetization in the determination 
of the spin accumulation and consider the steady solution. 
For the boundary condition of the spin accumulation, we 
artificially adopt 5m — at the top and bottom layers, 
ignoring any parasitic resistance. On the other boundaries, 
the natural boundary condition is employed. In the case, the 
simulated spin accumulation distribution is nonuniform and 
concentrates on the contact region. We adopt current density of 
0.01 mA/nm 2 where the top and bottom surface. The direction 
of current is aligned in the z-direction. 

Substituting the solution of Eqs. (f2|i and (01 into Eq. ((T), 
we evaluate the Hamiltonian of Eq. (Q~|i. The effective magnetic 
field for Si is given by dTL/dSi and the dynamics of S% is 
determined by the following Landau-Lifshitz-Gilbert equation, 



dt l 



7 



1 + a 2 



7 Si x 



dH 3 dH 

+ aSi x 
dS,. dSi 



(4) 



were 7 and a denote the gyromagnetic ratio and the Gilbert 
damping constant. The equation is numerically solved by a 
quaternion method[22|. The time step At and a are set at 
3.4 x 10~ 2 fs and 0.02, respectively. The number of sites, 
N s is about 10 4 . For simplicity, we use antiparallel boundary 
condition of the spins in order to simulate the experimental 
situation with a 180° domain wall. On the boundary, S's 
are fixed at (-1,0,0) on the top electrode and (+1,0,0) on the 
bottom electrode except for the boundary surface between the 
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Fig. 2. The relation between M and magnetic structures (a) and the time 
evolution of <j> for £ cx = 3nm(CoFe), 5nm(Permalloy), oo. The dashed and 
doted lines in the M x -M z and M x -M y planes denote the order parameters 
for the Neel and Bloch walls where M y and M z is zero, respectively. (The 
definition of M is exactly given by Eq.|5j. 



magnetic and insulator layers. The directions of spins, x, y 
and z are defined as shown in Fig. QI a )- 

To analyze the dynamics of the domain wall, especially 
the rotation of magnetic structures, as shown in Fig. |2ja) 
we calculate a vector, M which is the ferromagnetic order 
parameter scaled with its absolute value, 



M 



IM, 



op 



M, 



°p — TV X/ 1 ^' 



(5) 



(6) 



where N is the number of sites of local magnetizations. 
We show </> = arctan(M y /M z ) as a function of time t in 
Fig.|2jb). In the present setting the <f> almost linearly increases 
by spin transfer torque after enough time is spent. The angular 
velocities for all cases finally reach values of 1 xlO 2 Ghz. 
Because as shown in Fig. |2ja) the domain wall for <j> = mr 
and for <f> = 7r(l/2 + n), respectively, correspond to the Neel 
and Bloch walls, where n is any integer, the linear increasing 
of <j> suggests a magnetic structure rotation of the domain wall 
including the Neel and Bloch walls as discussed above. 

Figures [3ja)-(d) show snapshots of exchange energy density 
of the local magnetization for l cx = 3nm and oo with the 
Neel and Bloch walls. As shown in Figs. [3ja) and (b) for 
l cx =3nm the exchange energy density depends on whether the 
domain wall is the Neel and Bloch ones. In fact the maximum 
values of the exchange energy density of the Neel and Bloch 
are 0.08 and 0.11 [J /a 3 ], respectively. The difference of the 
exchange energy density in the magnetic structures of the 
domain wall suggests that the resistance depends on <fi and 
the resistance oscillation is induced by the rotational motion 
of (j> because both of the exchange energy and resistance due to 
potential drop of the spin accumulation is roughly proportional 
to squared inverse of the thickness of the domain wall[ 14|. On 
the other hand as shown in Figs. [3c) and (d) for l cx = oo 
the exchange energy density almost does not depend on a 
magnetic structure of the walls. In fact the maximum value of 
the exchange energy density is the almost same value of about 
0.1[ J/ a 3 ] for both the Neel and Bloch walls. In this case, the 
resistance can not oscillates during the rotational motion of (j>. 
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Fig. 3. The local magnetization configuration snapshots in the x-z plane 
including the center of the system and exchange energy density for (a) £ cx = 
3nm with a Neel wall, (b) Z ox = 3nm with a Bloch wall, (c) Z cx = oo with 
a Neel wall and (d) Z ox = oo with a Bloch wall. 



These results mean that the rotational motion driven by the 
spin-transfer torque induces a breathing mode which is an 
effective oscillation of the thickness of the domain wall. The 
breathing mode excitation induces the resistance oscillation of 
the system. As a result the microwave oscillation is expected 
because the oscillation of the resistance under a dc current 
corresponds an oscillation of voltage between the both ends 
of the contact. The breathing mode is excited only for the finite 
l cx because the effective thickness dependence on the magnetic 
structure of the domain wall is induced by the dipolar field. 
Thus we speculate that the magnitude of the voltage signal 
becomes small as increasing l ex and reaches zero at l cx = oo. 

Let us to directly examine the oscillation of the resistance. 
The electric field due to the spin accumulation is given by 



2D Q 



l + Pa- S(f) 5m{r) . (7) 



From the electric field we evaluate time dependence of 
the resistance. Figure. 0] shows the resistance ratio tmr = 
R(t)/Ro — 1 as a function of the time for l cx = 3, 5, oonm, 
here R(t) and Rq denote resistance and resistance without 
a domain wall. As speculated above the oscillation of tmr 
is observed for finite Z ex 's and the oscillation of tmr is 
absent for infinite l cx . The magnitude of the oscillation of the 
resistance increases with decreasing the exchange length. From 
the discussion so far we conclude that the enhancement of the 
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Fig. 4. The time evolution of rjyiR for l e 
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3nm(CoFe), 5nm(Permalloy), 



signal for short exchange length is because a strong dipolar 
field in a case of a short exchange length induces the large 
difference of the thickness of the domain wall between the 
Neel and Bloch walls. These results suggest that the magnitude 
of the microwave voltage signal is enhanced by adopting the 
material with short l cx if one design a domain wall spin-valve. 

Finally, we summarize our conclusions. We examined the 
dipolar field effect on the microwave generation in the spin- 
valve by solving simultaneously the Landau-Lifshitz-Gilbert 
and Zhang-Levy-Fert diffusion equations. Our results sug- 
gests that the dipole-dipole interaction crucially affects the 
microwave generation. By numerically analyzing the exchange 
length dependence, we showed that the efficient microwave 
generation needs the strong dipole-dipole interaction and the 
amplitude of the microwave voltage signal depends strongly 
on the exchange length. We also showed that the microwave 
voltage signal is induced by the breathing mode originating 
from the dipolar field. In order to develop the microwave 
device with the domain wall spin-valve with high efficiency, 
we suggest from the results that the materials with short 
exchange length are preferred so that the breathing mode with 
large amplitude is excited. 
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